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We examine structural and dynamical properties of quantum resonances associated with an 
avoided crossing and identify the parameter shifts where these properties attain extreme values. 
Thus the concept of avoided crossing resonance can be defined in different ways, which do not 
coincide in the general case. These definitions are described first at a general level, and then for 
a two-level system coupled to a harmonic oscillator, of the type commonly found in quantum op- 
tics. Finally the results obtained are exemplified and applied to optimize the fidelity and speed of 
quantum gates in trapped ions. 
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I. INTRODUCTION 

In this paper we study the definition of a resonance 
associated with an avoided crossing. Avoided crossings 
are quite common in many fields of nuclear, atomic, or 
molecular physics such as laser driven trapped ions [l|, Q , 
two level atoms coupled to a cavity mode or dia- 
magnetic hydrogen in magnetic fields R|- In the avoided 
crossing regions, two eigenvalues of the system approach 
first as a parameter ^ is varied, but then veer from each 
other. The "bare levels" of a zeroth order Hamiltonian 
do cross at a reference value but in the full Hamil- 
tonian a perturbation connecting them causes the split- 
ting. The eigenvalues also interchange their character 
along the avoided crossing: each of the two eigenvalues 
is dominated by different bare levels before and after the 
crossing. The resonance is also characterized by maxi- 
mal transition probabilities among the bare levels. As 
we shall see, however, the central loci of these phenom- 
ena do not generically coincide, so different "shifts" will 
be introduced. The shifts are not only due to the two 
bare levels directly implied, but to the "contamination" 
or influence of the rest of the levels. 

The article is organized in three increasing levels of 
concreteness: we shall first discuss formal general aspects 
in Sec. [ni then make a more specific analysis for Hamilto- 
nians commonly found in quantum optics which descibe 
a two-level system coupled to a harmonic oscillator in 
Sec. mil and finally exemplify and apply the results by 
optimizing the fidelity and speed of quantum gates for 
trapped ions in Sec. IIVI The paper ends with Sec. |Vl a 
brief discussion of the results obtained, and a technical 
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Appendix. 

II. GENERIC AND FORMAL ASPECTS 

Consider the general Hamiltonian 

H = HoiO + V{0, (1) 

where V{^) is a small perturbation of the non-perturbed 
(bare) Hamiltonian Hq. Both parts may depend on an 
external parameter ^ and the eigenenergies and eigen- 
states of Hq (bare energies and states) are supposed to 
be known, 

i/o|a) =e„|a). (2) 

with = huja ■ Assume also that the energy levels corre- 
sponding to two given bare states \a) and \b) cross each 
other at {^o,Eo) in the {^,E) plane (Fig. [T^, dashed 
lines). This defines the bare or unperturbed resonance 
at ^ = ^0- If both states are connected by the perturba- 
tion, Vab — (a\V\b) ^ 0, the crossing between these levels 
will become an avoided-crossing. 

A. Structural Shift 

Let us now assume that the two energy levels and 
£(, are close to each other but far from other levels, i. e., 
a well isolated avoided crossing. The perturbed energy 
levels may be exactly given by an effective Hamiltonian 
[IjUjQ in the subspace spanned by the states \a) and \b) 
given by 
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FIG. 1: (a) Schematic plot of the energy levels around the 
resonance between the \a) and bare states. Non-perturbed 
(dashed lines) energy levels ta and eb cross each other at the 
bare resonance ^ = ^o. The perturbation V shifts the position 
of this resonance (structural definition) to ^ = (solid lines), 
(b) According to a dynamical definition the resonance is in- 
stead at 5 = ^_D, where the state-flip probability is maximum. 
The dynamical shift A_d is the separation between both defl- 
nitions, i. e., Ad = — £,s- (c) The change of character of 
the dressed energy levels occurs at ^n. 



where the effective detuning 5 is defined as 
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5 = - (Wfc - -I- nb - raa) 



(4) 



and Tap are the matrix elements of the level shift operator 
divided by fi,, 



r{E) 



R{E) 



R{E) = PVP + ^PV 



Q 



E-Hq 
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V] P, (6) 



with P = \a){a\ + \b){b\ and Q = 1 - P. Even though it 
is 2-dimensional, the effective Hamiltonian ^ contains 
information about the whole Hilbert space via the op- 
erator Q, the projector onto the non-resonant subspace. 
iJe/ / is an implicit Hamiltonian, since it depends on the 
exact perturbed energy E through R{E). A first approx- 
imation is to evaluate Hej / at the nonperturbed position 
of the resonance, i. e., at ^ = and E — Eq. Further 
corrections can be obtained by iteration [3]. 

The position of the resonance may be defined to be 
at the center of the anti-crossing 0, [^, which may be 
found at the extrema of the dressed energy levels, i. e., 

C = ?S in Fig. [T^. This is also the point of minimal 
splitting. Using the expression of the eigenergies of Heff, 



|r.6l')"' = 0, 



or more explicitly 



'4 



d\rab\ 



= 0, 



(7) 



(8) 



see Fig. [T^. 

The discussion so far refers to an energy level structure 
criterion so that As = — is a "structural shift". A 
different definition, based on the dynamics of the system, 
will be discussed in the following. 

Before, note that if PVP does not contribute, A5 is 
a generalization of the Bloch-Siegert (BS) shift, first de- 
scribed by Bloch and Siegert while studying the failure 
of the Rotating Wave Approximation [3, 7], and due en- 
tirely to the effect of non-resonant transitions involving 
levels in the subspace Q. 



B. Dynamical Shift 

Let us now prepare the system in the \a) bare state and 
look for the probability to find it in |6). We could define 
the resonance in a dynamical way as the value of ^ where 
the state-flip probability is maximal. For the dynamics 
governed by the effective Hamiltonian the state-flip 
probability is given by 
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where fl — (S'^ + |fa6p)^^^- If a 7r-pulse (defined by Qt 
tt) is applied, this probability shows a maximum for 6 
0, at ^ = defined by the condition 



S = LJb - UJa + Tbb - Taa = 0, 



(10) 



see Fig. [TJd. 

The dynamical definition of the resonance may be also 
understood in terms of the change of the character of 
each dressed energy level. Consider the eigenstates of 

Heff 



\e+) = cos^e'^^^\a)+sm-e~'^^^\b), 

= -siJ-e^^/^\a)+cos^-e-^*/^\b), (11) 



^ The exact dynamics in the P-subspace is given by a non- 
Markovian equation with a memory kernel, and H^f / provides 
the Markovian approximation. 



3 



where tanO 



\rab\/S and rab 



Fab 



e*"^. The "char- \a,na) and \b,ni,) cross each other at 



acter" of the dressed states around resonance is given 
by Pa.,± = |(a|e±)P, the projection of the bare state 
= a,b) onto the dressed state |e-|-). Then we have 

that 



Pa+ = Pb- = COS -, 



Pa- = Pb+ = sm 



(12) 



The character change of the dressed levels is centered at 
9 = 7r/2, when both projections are equal, Pa+ = Pa-- 
This corresponds to (5 = 0, and coincides with the dynam- 
ical resonance condition (fTO|) . see Fig. [TJ:. may thus be 
defined independently of the 7r-pulse condition. Numer- 
ical calculations diagonalizing full n-dimensional Hamil- 
tonians (with n large enough to assure convergence) con- 
firm this result, and examples are provided in Sec. IIVI 
below. 

To summarize the main results obtained so far: the 
resonance location may be defined in different ways: 
a "structural" criterion, Eq. ([8]), gives ^s, where the 
splitting between the dressed energy levels is minimal, 
whereas the "dynamical" criterion gives ^_d, defined by 
the condition ^TU^ . where the state-flip probability is 
maximal, and do not coincide in the general case, 
since the two conditions cannot be simultaneosly full- 
filled as long as |rab| depends on The two values are 
separated by the "dynamical shift" A/j = — ^5, pro- 
portional to d\rab\'^/d^. 



^0 = (nb - na)uj, 

huj 

Eo = [ua + nb) — , 



(15) 
(16) 



see Fig. [T^, dashed lines. As pointed out in the previous 
section, the perturbation will not only split the energy 
levels but it will also shift the position of the anti-crossing 
from ^0 to ^S: see Fig. [T^, solid hues. It may also happen 
that both levels are shifted {raa, fbb 7^ 0) but not splitted 
(Tab = 0), and thus the crossing remains permitted, as in 
Fig. [1^, dotted lines. The perturbed energy levels will 
be described by the effective Hamiltonian with the 
effective detuning (|4]), 



(^0 - C + nb - Taa) ■ 



(17) 



The position of the resonance according to structural and 
dynamical criteria will be determined by the conditions 

m and (nni), 



is = Co + ^-fcb - -f Ad(^), 

S,D = Co + rbb - raa, 



Ad(0 
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(18) 
(19) 

(20) 



where, as a first approximation, all the matrix elements 
of r and their derivatives are evaluated at Co- 



lli. TWO-LEVEL SYSTEM COUPLED TO A 
HARMONIC OSCILLATOR 

One of the simplest cases in which the difference be- 
tween the structural and dynamical definitions of a reso- 
nance can be observed is a two level system coupled to a 
harmonic oscillator and driven by some external field, as 
frequently found in quantum optics. The Hamiltonian, 
in an appropriate interaction picture takes the form 



H 



noja'a + —a^ 



v^(e), 



(13) 



where V^(C) is assumed to be a small perturbation of the 
unperturbed Hamiltonian iJp = Titjata + ?iCfz/2- and 
a are the usual creation and annihilation operators for 
the harmonic oscillator, and the Pauli atomic inversion 
operator is = |a)(a| — V may also depend on 

Oi with i = x,y, z. The bare energy levels of this system 
are 



nhuj 



nhuj 



2 ' 



(14) 



n being the vibrational quantum number n = 0, 1, 2 . . .. 
The energy levels corresponding to two given bare states 



IV. TUNING QUANTUM GATES FOR 
MAXIMAL SPEED 



The efficient physical implementation of quantum 
gates and quantum information processing is a major 
goal for different fields of physics. The approaches based 
on ions in a linear trap pioneered by Cirac and ZoUer Q 
are among the most developed, and have become a work- 
ing horse to test basic quantum information processing 
0. In addition, the formalism is very similar or even 
equal in some limits to the one applied in other systems, 
such as cavity QED [l3|, or superconducting qubits [ll|. 

Quantum gates based on trapped ions illuminated by 
lasers can be speeded up considerably by adjusting the 
laser to the exact position of the resonance [ij] ■ In these 
quantum gates one is interested in obtaining the maxi- 
mum fidelity, that is, the maximum transition probability 
from one state to another. The dynamical shift defined 
above plays then an important role, since the laser pa- 
rameters have to be adjusted to Cd and not to ^5, unless 
of course they coincide. 

Let us now consider an effectively ID trapped ion inter- 
acting with a classical field in a laser adapted interaction 
picture and after applying the optical Rotating Wave Ap- 
proximation (RWA). It is described by the Hamiltonian 



hA h^R 

H 

2^2 



H.c 



(21) 
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where A = lol — wo is the detuning (laser frequency 
minus transition frequency between levels \g) and |e)), 
Tj = (wa/wt)^^^ is the Lamb-Dicke (LD) parameter, and 
fz = |e)(e| — 1 5) (5 1; is the recoil frequency of the ion 
and ^Ir is assumed real without loss of generality. 

A. Stark Shift gate 

The dynamical shift is clearly observed in the so-called 
"Stark shift gate" , proposed by Jonathan, Plenio and 
Knight [3] ■ Looking for faster quantum gates, they pro- 
posed a scheme where high intensity lasers overcome the 
slowness problem of the Cirac-Zoller (CZ) Q and re- 
lated gates, where low intensity lasers (in the sense that 
rifl ^ ut) limit the gate velocity. Setting the laser fre- 
quency on resonance (A = 0) and the laser intensity so 
that the Rabi frequency and the trapping frequency 
ujt coincide (first "Rabi Resonance" — ^t), the 

splitting of the dressed states coincides with one vibra- 
tional quantum, so we can expect avoided level crossings 
and fast and efficient vibronic transitions. This gate is 
thus based on a double resonance condition for both laser 
frequency and intensity. 

The Hamiltonian given in Eq. (j2ip is written in the 
basis of the eigenstates of CTz {|5),|e)} which form the 
computational basis in the CZ gate scheme. The Stark- 
Shift gate, on the other hand, works in the {|±)} basis, 
where |±) = (I5) ± |e))/\/2 are the eigenstates of ax- It 
is therefore convenient to write H in this new basis, 

t 

+ — Y~ [cos 77(0 -I- a')az + i sin7y(a -I- a^)((7+ — (7_)] , 

where we have redefined the Pauli operators in the {|±)} 
basis according to CT+ = |+)(— | =d'^_,d'x = a++a-,az — 
|+)(+| — |— )(— |. These tilde operators are related to the 
usual Pauli operator by the transformation 

<yz = -^x, (22) 

(T+ = - (tr^ + (7-+ - = (t1, (23) 
(7x = a^. (24) 

Both Hamiltonians H and H are completely equivalent, 
only a change of basis has been applied. In order to 
allow for lasers with arbitrary intensity we consider the 
LD parameter to be the perturbative parameter. Then, 
we may split H as Hq + V{ri) with 

Ho = hujTa'a^ 2~^'^' 

^(.v) = — [(cosci — 1) (T2 + i sind ((T+ — (T_)](25) 

where a = r]{a -I- a^) and where A = as required by 
the SS-gate (the perturbation vanishes in the LD limit, 




FIG. 2: (Color online) Stark Shift gate scheme: (a) Energy 
levels of a trapped ion coupled to a resonant (A = 0) laser 
around the first Rabi Resonance where the Stark Shift gate 
transition |-|-,0) | — ,1) takes place. The non perturbed or 
bare (r; = 0.0) energy levels e+,o and e-,i cross each other at 
(yiE)o, dashed lines. The perturbed energy levels (77 = 0.3) 
form an avoided crossing centered at {Q.r)s, and shifted from 
{yiRjn by As, solid lines. The dash-dotted (red) lines repre- 
sent the derivative of the perturbed energy levels with respect 
to f2_R. They cross each other at {Q,ii)s, where the deriva- 
tive vanishes (zero slope), (b) When a vr-pulse (defined as 
•qO.Rt — tt) is applied, the state-flip probability from |-|-,0) 
to | — , 1) shows a maximum at {flR)D, shifted from {flR)s 
by the dynamical shift Ad. (c) The change of character of 
the dressed energy levels occurs at {i}[i)D- AH calculations 
have been done by numerically diagonalizing the full Hamil- 
tonian (|2ip to all orders in rj and including a large number of 
vibrational states (n = 25). 



V(rj = 0) =0). This SS-gate Hamiltonian has the same 
form of the general Hamiltonian p3p . so the formalism 
described in Sec. IIIII may be applied by redefining \a) = 
\ + ,n), \b) ^ h,n+l> and ^^flR. 

B. The stark-shift gate error and velocity 

Let us now consider the Stark-Shift gate interaction 
between computational states and | — ,n-|-l). If the 

system is prepared initially in the state |-f,?^), an ideal 
Stark-Shift operation would flip the state of the system 
toj — 1). We define the gate error as the imprecision 
[12| given by 

e«= [l-P|+,„>^|-.„+i)]'^'. (26) 

Numerical simulations with a large number of vibrational 
states show that this error is a linear function of the LD 
parameter 77. To achieve the smallest error, one needs to 
correct the resonance position. We may now apply the 
formalism described in Section Hill to the Stark-Shift gate 
Hamiltonian (pS)) . with ^ = O^^, see Appendix [XI to find 
that the positions of the structural and dynamical reso- 
nances for the \+,n) | — ,n-|-l) transition are shifted 



FIG. 3: (Color online) Error of the quantum gate after a tt- 
pulse as the perturbation becomes stronger (a) in the SS-gate 
scheme and (b) in the CZ scheme with = 0.1. Solid (black) 
lines: the laser is tuned to the bare resonance without any 
correction. Dotted (red): the laser is tuned to the structural 
resonance. Dashed (blue): tuning the laser to the structural 
resonance in the SS-gate gives a larger error than the uncor- 
rected gate, this shows the importance of the "dynamical" 
shift. In the CZ scheme (b) the dynamical shift is negligible 
so the dotted (red) and dashed (blue) lines coincide. In both 
gates, the dynamics with the effective Hamiltonian would give 
a zero error at the dynamical resonance. The discrepancy 
with the exact calculation is due to the non-Markovian char- 
acter of the true evolution. 



FIG. 4: Cirac-ZoUer interaction (a) First few energy levels of 
a trapped ion at the first blue sideband, where the states !<;, 0) 
and |e, 1) are on resonance. Point A marks the bare nominal 
resonance at A = lot when the laser is turned off (f2_R = 
0), while points Bl and 52 define the structural resonance 
(the position of the anti-crossing) if the laser is turned on 
(yia/iOT ~ 0.3). (b) State-flip probability from \g,0) to |e, 1) 
after a vr-pulse (defined by r^Qnt = tt) has been applied; the 
maximum at C defines the resonance location dynamically. 
In the CZ interaction both structural dynamical definitions 
coincide, there is no dynamical shift, i. e.. Ad = 0. All 
calculations have been done by numerically diagonalizing the 
full Hamiltonian (|2H) to all orders in 77 and including a large 
number of vibrational states (n — 20). 



from the bare position fi^^ = cut according to 

1 

—J 

3 



ujt — -r]'^ujT{n + 1), 



(27) 
(28) 



Tuning the Rabi frequency to {^r)d does not change the 
dependance of the error with 77, which remains linear, but 
improves by a factor of 3 the obtained error, see Fig. 
[3^, while tuning to (l^fl)s the error is worse than for the 
bare resonance. 

Since the Rabi frequency, and thus the gate velocity, 
are limited by the gate error the threshold value 
et of the imprecision of the gate will limit the processor 
speed. For a given required gate precision, the velocity 
of the SS gate (expressed as the inverse of the 7r-pulse 
duration) will be limited by the condition 



1 



(29) 



This bound can be improved by a factor of ~ 3 by opti- 
mization of the intensity as discussed above. The speed 
bound is better by I/77 than the bound of the gate of 
Monroe et al. [11] or the corrected CZ gate p^ . 



C. Comparison with the CZ gate 

It is interesting to compare the SS scheme to the CZ 
gate scheme, where the perturbative parameter is no 
longer the LD parameter. In this case, the small param- 
eter is the Rabi frequency (low intensity lasers), "small" 
meaning <C cot- Moreover, the computational basis 
in this scheme is the bare \g), |e) basis (for low intensity 
lasers the dressed states may be approximated by the 
bare states). The Hamiltonian (PT|) may then be parti- 
tioned as 



Hq — hujTCi'a —CTz 



~2 



V 



H.c 



(30) 



Applying once again the effective Hamiltonian formalism 
described in Section Hill (redefine \a) = \g), \b) = |e), and 
^ = A), it is found that both structural and dynamical 
resonance definitions coincide, see Fig. H) Thus, the CZ 
interaction shows no dynamical shift, Ajj — 0. In order 
to speed up the CZ gate, correcting for the structural 
shift is enough see Fig. ^jp. Explicit expressions for 
the structural shift for a CZ type of interaction (the so- 
called Vibrational BS shift) have been provided elsewhere 



6 



V. DISCUSSION 

Two different definitions of a quantum resonance as- 
sociated with an avoided crossing have been provided, 
which do not coincide in the general case. As an ap- 
plication, we have shown that, for the same precision 
requirements, the speed of the SS gate is of order l/rj 
times larger than the speed of the corrected CZ gate, 
but its optimization implies the new concept of tuning 
the laser intensity to the combined effect of dynamical 
and structural shifts. We expect that this finding will 
have repercussions in other gates affected by light shifts 
0, [3 J and in related physical systems (cavity QED, 
superconducting qubits in which a similar Hamil- 

tonian structure and fast high intensity transitions are 
considered [l^]- The distinction between dynamical and 
structural shifts will also be relevant for metrological ap- 
plications fT7| , as in atomic clocks and other interferom- 
eters. 
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APPENDIX A: STRUCTURAL AND 
DYNAMICAL SHIFTS FOR THE SS GATE 
INTERACTION 

The matrix elements of the level shift operator R and 
thus the structural or vibrational BS shift for the CZ type 
of interaction (low intensity lasers) are provided in 
In order to obtain explicit expressions for the dynamical 
and structural resonances defined in Sec. [Till for the SS 
gate Hamiltonian, we calculate the matrix elements of 
the level shift operator using the SS-gate Hamiltonian 
given in Eq. ((25)) . 



2 (^»+."+ 



^2 



E 



_|r|-,n_ 



E 



E 



r+^{E,nn) - {+,n+\r\- 



+ E 

SIr 
* 2 



E 



(Al) 



with Cn.n' — cos a|rt') = Rc(fin_n') and Sn,n' = 

U,ji{n\ sindirt') = Im {^n.n') being the real and imaginary 
Bxt of the coupling strenghts ^n,n' — {n\e"'^°'^°'''^\n') 
and where the non-perturbd energy levels are given 



[1 

by 



e± n = nh(jjT ± 



(A2) 



As pointed out in 0, Rabi resonances will be well de- 
fined when the corresponding avoided crossings are well- 
isolated, i. e., when the ion is confined within the LD 
regime 77 ^ 1 [2] . In that case the matrix elements of R 
to leading order in t] will only have contributions from 
the "nearest neighbour" vibrational states. The matrix 
elements are thus simplified to 



= r,'nnT±{E,nR), 



'+- 



«?y— V"'-On+,n_-l, 



(A3) 
(A4) 



with 



1 



2 

4 V£^ 



2 

n± 



n± 



ET,n±-l 



E- 



(A5) 



and where any term with a zero in the denominator has to 
be omitted since it is excluded from the sum. Since these 
matrix elements are part of the implicit Hamiltonian ([3|) , 
the first approximation is evaluated at the non-perturbed 
position of the resonance, i. e., at -E = -Eq = ("■+ + 
n-)hu!T/2 and Hr^ = fJ^^ = (n_ — n-f-)ujT- 

With these expressions it is easy to find the structural 

shift for an arbitrary resonance, A5 = r — r++ + Ad, 

see Eq. ([T5|) . In particular, we find for the first few 
resonances 



(A5 


n,n+l ^ 


1 , 

—^{n + 


1) V^UJT, 


(A5 


n,n+2 — 


-i(2n- 
3 ^ 


h 3) 


(A5 


n,n+3 


-g(n + 


2) V^tOT, 


(As 


n,n+4 — 


2 

-— (2n 
15 ^ 


+ 5) 77^ wt 



(A6) 



For the particular case of the SS-gate transition \+,n) ^ 
| — , n-l- 1) discussed in Sec. IIVI we have, to leading order 
in ?7, 



r++(ri_R = lot) 



|?7V(3n-|-2), 



1 



r — (flR^ujT) = -^V cjT(3n-f4), 
o 

|2 



d\r+- 



R 



^ ^^"^(n + l). (A7) 



The positions of the resonance according to structural 
and dynamical definitions for this SS-gate transition (as 
well as the value of the dynamical shift) are now readily 



7 



calculated using Eqs. ((T8l) and (fTQ]) . which corresponds to a dynamcial shift | A^ij = rj^u}T{n + 

1). Note that the derivatives of the diagonal terms r__ 
9 1?'^ 1^ and r__ for computing the dynamical shift from Eq. ([^0]) 



(i^fl)s i^/j + T r++ 2 ^j-^^ have been neglected, since they only contribute with 77* 

I terms. 
= - 7'7^WT(n + 1), (A8) 



WT + ^'7^WT(n+ 1), (A9) 
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